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ABSTRACT 

The  distortion  of  a  two-dimensional  bubble  (or  drop)  in  a  corner  of  angle 
6,  due  to  the  flow  of  an  inviscid  incompressible  fluid  around  it,  is 
examined  theoretically.  The  flow  and  the  bubble  shape  are  determined  as 
functions  of  the  angle  <5,  the  contact  angle  0  and  the  cavitation  number 
Y.  The  problem  is  formulated  as  an  integrodif ferential  equation  for  the 
bubble  surface.  This  equation  generalizes  the  integrodif ferential  equations 
derived  by  Vanden-Broeck  and  Keller1'2,  The  shape  of  the  bubble  is  found 
approximately  by  using  the  slender  body  theory  for  bubbles  presented  by 
Vanden-Broeck  and  Keller2,  When  Y  reaches  a  critical  value  Yq(0,<5), 
opposite  sides  of  the  bubble  touch  each  other.  Two  different  families  of 
solution  for  Y  <  YQ  are  obtained.  In  the  first  family  opposite  sides  touch 
at  one  point.  In  the  second  family  contact  is  allowed  along  a  segment.  The 

methods  used  to  calculate  these  two  families  are  similar  to  the  ones  used  by 

3  4 

Vanden-Broeck  and  Keller  and  Vanden-Broeck  . 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  mixing  of  two  fluids,  a  drop  or  bubble  of  one  fluid  will  be 
distorted  and  possibly  split  into  smaller  parts  because  of  the  flow  of  the 
other  fluid  around  it.  In  order  to  study  this  phenomenon,  we  consider  a  two- 
dimensional  bubble  (or  drop)  in  a  corner  of  angle  6  (see  Figure  1).  The 
bubble  is  distorted  by  the  flow  of  an  inviscid  incompressible  fluid  around 
it.  We  take  into  account  the  surface  tension  o  at  the  interface,  but  we 
ignore  the  flow  inside  the  bubble,  assuming  the  pressure  is  a  constant  pb 
throughout  it.  The  flow  and  the  bubble  shape  are  determined  as  functions  of 
the  angle  <5,  the  contact  angle  3  and  the  cavitation  number  Y. 

The  problem  is  formulated  as  an  integrodifferential  equation  for  the 
bubble  surface.  This  equation  can  be  solved  numerically  by  using  the  scheme 
derived  by  Vanden-Broeck  and  Keller  9  •  In  addition  the  shape  of  the  bubble 
can  be  found  approximately  by  using  a  slender  body  theory. 

For  Y  =*  w  and  3  =  Figure  1  represents  half  a  free  bubble  in  an 
uniform  stream.  For  Y  «  8  3  p  Figure  1  represents  a  quarter  of  a  free 
bubble  in  a  straining  flow.  Explicit  results  are  presented  in  these  two 
important  cases. 

When  Y  reaches  a  critical  value  YQ(B,<5)  opposite  sides  of  the  bubble 
touch  each  other.  Two  different  families  of  solution  for  Y  <  YQ  are 
obtained.  In  the  first  family  opposite  sides  touch  at  one  point.  In  the 
second  family  contact  is  allowed  along  a  segment. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


BUBBLE  IN  A  CORNER  FLOW 


Jean-Marc  Vanden-Broeck 

1.  Introduction  and  formulation 

We  consider  the  steady  potential  flow  aroand  a  gas  babble  or  liquid  drop  in  a  corner 
of  angle  6.  The  contact  angle  is  denoted  by  0  (see  Figure  1)*  We  shall  write  "babble" 
to  mean  either  bubble  or  drop.  We  take  into  account  the  surface  tension  o  at  the 
interface,  but  we  ignore  the  flow  inside  the  bubble,  assuming  that  the  pressure  is  a 
constant  throughout  it. 


Figure  1.  Sketch  of  the  bubble  and  the  coordinates 

In  order  to  formulate  this  problem  we  assume  that  the  complex  potential  without  the 
6a  n  /6 

bubble  is  —  (x  +  iy)  ,  where  a  is  a  constant  and  x  and  y  are  Cartesian 
*  '*>ordi  nates. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  .  This 
material  is  base  1  upon  work  supported  by  the  National  Science  Foundation  nnd-»r 
Grant  No.  MCS -7427062,  Mod.  1. 


as  the  unit  length  and 


r  2C\ 

We  Introduce  dimensionless  varia1  es  by  choosing  [ — -J 

pa 


6 


as  the  unit  velocity.  We  also  introduce  the  dimensionless  potential  b$  and 


u-6 

r  20  ^ 
o(— ) 

po 

stream  function  Here,  b  >  0  is  a  dimensionless  constant  to  be  chosen  so  that 

$  -  and  t  -  -  j  at  the  stagnation  points  on  the  walls  y  ■  0  and  y  -  xtan6, 
respectively.  We  denote  the  streamline  along  the  two  yalls  and  along  the  bubble  boundary 
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by  ♦  -  0.  In  these  variables  b($  +  !♦)  ~  ~  (x  +  iy)  at  Infinity  or,  equivalently 


♦  iy  ~  Onr)  (♦  + 


O) 


at  infinity. 

The  flow  occupies  the  region  ^  >  0  of  the  K  t  plane,  and  the  bubble  boundary 
corresponds  to  the  segment  -  <  $  <  -j  of  the  axis  ^  *  0.  The  problem  of  finding  the 

flow  consists  of  determining  x  +  iy  as  an  analytic  function  of  ♦  +  ity  in  the  half  plane 
>  0  satisfying  Equation  (1)  at  infinity.  Then  the  bubble  surface  is  given  by  setting 
t  =  0  in  x(t  +  it)  and  y(<t  +  it)  and  letting  t  range  from  -  to  The  contact 

angle  conditions  require  that  the  bubble  surface  meets  the  walls  at  the  angle  0,  which 
yields 


tan  0  as 


t 


2 


(  -tan(0  -  6)  as  $  ♦  -  . 

On  the  bubble  surface  the  pressure  in  the  fluid,  which  is  given  by  the  Bernoulli 
equation,  must  differ  from  by  crk,  where  k  is  the  curvature  of  the  interface, 
leads  to  the  boundary  condition 


(2) 


This 


pg  -  «  pb  -  <?k  on  -  J  M  t  -  0  .  (3) 

Here,  pg,  p  and  q  are,  respectively,  the  stagnation  pressure,  the  density  and  the  speed 
of  the  fluid  outside  the  bubble.  In  dimensionless  variables  (3)  becomes 


-2- 


♦  -  0  , 


(4) 


q2  »  k  -  Y  on  ~  J  <  +  <  5' 
where  Y  is  the  cavitation  number  defined  by 


Y 


(5) 


The  problem  can  be  further  simplified  by  requiring  the  bubble  to  be  symmetric  about 

Y 

the  line  y  *  x  tan  —  •  This  implies  that 

y^ ($,°)  *  y^(-^,0)cos6  -  x^(-^,0)sin6,  0  <  $  <  j  .  (6) 

By  using  Equation  (6)  we  can  restrict  our  analysis  to  the  interval  0  <  ♦  < 


-3- 


«nw-  <>k-  i. 


2m  Reformulation  as  an  integrodlf ferential  equation 


It  is  convenient  to  reformulate  the  boundary  value  problem  as  an  integrodlf ferential 
equation  by  considering  the  function 


1- 

(9  +  i^) 


which  is  analytic  in  the  half  plane  9  >  0  and  vanishes  at  infinity  as  a  consequence  of 
Equation  OK  Therefore#  on  9  ■  0,  its  real  part  is  the  Hilbert  transform  of  its 
imaginary  part.  The  imaginary  part  vanishes  on  ^  »  0,  | $ |  >  and  therefore  the  Hilbert 
transform  yields 


1- 

♦ 


6  1 
V(^o»  -  (Vf  ■ 


1  V2  (♦') 
i r  * 


♦'  -  ♦ 


<*♦' 


♦*/ 


(-♦ 


>'mh 


-1/2 


■y^(<fr* #0)cos6  +  x^($',0)ain6] 


d$ ' 


(7) 


We  now  use  the  symmetry  condition  (6)  to  rewrite  (7)  in  the  form 


x^(*,0) 


_6 

ir 


-1 


i-i 


(8) 


Next  we  express  the  boundary  condition  (4)  in  terms  of  x.  and  y.  noting  that 

9  9 

2  2  2  2  -1 

q  -  b  (xA  ♦  y.)  •  Then  (4)  becomes 

9  9 


2  2 

*♦ +  y* 


yit 


-  «£li_ 


<*;  *  v2s/2 


Y, 


1*1  <  J,  * 


(9) 


Now  (8)  and  (9)  together  constitute  a  nonlinear  integrodif ferential  equation  for 

y^($)  in  the  interval  0  <  ♦  <  -j#  ^  *  0.  The  contact  angle  conditions  (2)  complete  the 

formulation  of  the  problem  for  yA (♦,(>)  and  b. 
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For  y  =  n ,  the  equation  defined  by  (8)  and  (9)  reduces  to  the  integrodif ferential 
equation  derived  by  Vanden-Broeck  and  Keller1,  The  particular  case  8  -  j  represents  half 
of  a  free  bubble. 

For  y  =  8  *  the  equations  (8)  and  (9)  reduce  to  the  integrodif ferential  equation 
derived  by  Vanden-Broeck  and  Keller^.  This  case  represents  a  quarter  of  a  free  bubble  in  a 
straining  flow. 

The  integrodif ferential  equation  defined  by  (8)  and  (9)  can  be  solved  numerically  for 
arbitrary  values  of  8#  y  and  6  by  using  the  numerical  procedures  described  by  Vanden- 
Broeck  and  Keller1'^. 

In  the  next  section,,  we  shall  find  the  shape  of  the  bubble  approximately  by  using  the 
slender  body  theory  for  bubbles  presented  by  Vanden-Broeck  and  Keller1. 
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3.  Slender  body  approximation 

For  Y  large  the  bubble  bends  to  an  arc  of  a  circle  of  radius  Y  1 •  As  Y  decreases 

numerical  solutions  show  that  the  bubble  elongates  in  the  direction  of  the  line  which 

bisects  the  angle  between  the  two  walls.  Then  it  develops  a  horn  or  spike  which  large 

curvature  near  its  end.  Finally  when  Y  reaches  a  critical  value  Yq(0*6),  opposite 

sides  of  the  bubble  touch  each  other.  For  0  <  the  contact  point  is  at  x  *  y  =  0.  For 

0  >  -j  the  contact  point  is  away  from  x  =  y  *  0.  Typical  profiles  for  6  *  x  and 

6  *  0  «  can  be  found  in  Vanden-Broeck  and  Keller1'2.  These  profiles  were  obtained  by 

solving  numerically  the  integrodif ferential  equation  of  Section  2. 

For  Y  ~  Yo(0,«>  the  bubble  is  slender.  Therefore  we  shall  use  the  slender  body 

theory  for  bubbles  presented  by  Vanden-Broeck  and  Keller1  to  get  an  approximate  description 

of  the  flow  around  the  bubble.  In  the  lowest  order,  the  flow  about  a  symmetric  slender 

bubble  is  approximated  by  the  flow  about  a  rigid  plate  lying  along  the  center  line  of  the 

bubble.  In  the  present  case  the  center  line  of  the  bubble  consists  of  a  straight  segment 

of  some  length  a  lying  along  the  line  y  *  x  tan6/2.  We  introduce  the  coordinates 

x',y'  (see  Figure  1)  and  find  the  potential  b$(x',y')  of  the  flow  about  these  plates 

x 

6  1 

requiring  that  at  infinity  b($  +  i\Jj)  **  —  (x  +  iy)  .  Evaluating  the  potential  on  the 

IT 

plate  y'  ~  0,  x'  >  0  we  obtain 


2x  2x  1/2 

6  r  T  "IT-. 

b$(x',0)  -  —  [a  -  x'  ] 

By  differentiating  (10)  we  find  that  the  flow  speed  q  on  the  plate  is 


(10) 


2x  2x  2x  -1/2 

s  —k  -j\ 

q(x',0)  -  x'  (a  -  x‘  )  ,  x'  >  0  .  (11) 

Before  using  q  to  get  the  bubble  shape,  we  shall  determine  the  length  a.  We  do  so  by 

requiring  the  suction  force  F,  exerted  by  the  flow  on  the  end  of  the  spike,  to  balance 

2 

the  surface  tension  2a.  As  we  see  in  Ref.  5  [p.  412,  Eq  (6.5.4)],  F  *  xpA  /4.  Here 
A,  is  the  coefficient  in  the  expansion  b$  Ar1^2  cos  ~  in  terms  of  polar  coordinates 
with  their  origin  at  the  end  of  the  plate.  Upon  setting  F  =  2a  and  introducing 


dimensionless  variables  we  obtain 


(12) 


2n 


From  (10)  we  find  A2  «  —*  a 

It 


-1 


A2  =4  • 


so  (12)  yields 


(!) 


_ 6 

2* -6 


(13) 


We  next  use  (11)  for  q  in  (4)  and  approximate  the  curvature  k  by  -n^.^tx')* 
Here  the  equation  of  the  bubble  is  y*  *  n(x*).  Then  (4)  becomes 


n  ,  ,  »  -x' 
x'x' 


4ff 

“T 


2ir 


2*  -1 


(a  °  -  x.  °) 

At  the  end  of  the  spike  we  require 

n(a)  =  0 

In  addition  the  contact  angle  condition  yields 


n* (a)  *  -tan ( B  -  — ) 


Here  a  is  defined  by  the  equation 


<$ 

n(a)  *  atan  —  • 


(14) 


(15) 


(16) 


(17) 


The  function  n(x')  is  easily  obtained  by  integrating  (14)  twice  with  the  auxiliary 
conditions  (15)  and  (16).  In  the  particular  case  6  =  n,  the  result  of  the  integration  is 
n(x*)  *  (a2  -  x' 2 ) (Y  -  1  )/2  -ya(a  +  x’  )log(a  +  x'  ) 


-  -j  a(a  -  x')log(a  -  x*)  +  a2log2a  +  (a  -  x')tan(B  -  y)  .  (18) 

For  6  a  w,  (13)  becomes 

a  =  j  .  (19) 

Vanden-Broeck  and  Keller1  have  shown  that  the  approximate  solution  (17),  (18)  is  in  fair 
agreement  with  the  exact  numerical  solution  of  (8)  and  (9)  for  3  ~  ir  and  Y  **  Y^tB,*). 

For  Y  <  YQ(8,6)  ( 1 4) -( 1 6)  yield  unphysical  profiles  in  which  opposite  sides  of  the 

bubble  cross  over.  In  the  next  two  sections  we  construct  physically  acceptable  families  of 
solutions  for  Y  <  YQ(3/6 ).  we  shall  present  these  results  in  the  important  particular 
case  6  *  it  . 
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Solution  with  one 


»int  of  contact 


J22 

To  obtain  solutions  for  Y  <  w ^ 3 )  we  require  the  free  surface  to  be  in  contact 
with  itself  at  one  point.  Then  the  bubble  contains  a  small  sub-bubble  near  its  tip  (see 
Figure  2).  We  denote  by  c  the  x*  coordinate  of  the  contact  point. 

We  describe  the  profile  of  the  bubble  by  the  equations  y'  *  h^(x')  0  <  x'  <  c 

and  y*  =  n^fx')  c  <  x*  <  a.  Then  by  symmetry  we  have 

Vc>  *  n2(c)  *  0  {2{ 

nj(c)  *  n^(c)  *  0  .  (2 


X* 


Figure  2.  Profile  of  the  bubble  with  one  point  of  contact  for  Y  *  YQ  ~  -1,7  and 
0  ■  2*/3.  The  vertical  scale  is  the  same  as  the  horizontal  scale.  The 
cavitation  number  in  the  sub-bubble  is  equal  to  Y  . 


Figure  3.  Profile  of  the  bubble  with  one  point  of  contact  for  Y  =  -4.5 

and  B  =  2*/3.  The  vertical  scale  is  the  same  as  the  horizontal 
scale.  The  cavitation  number  in  the  sub-bubble  is  y  -  -0.6. 

The  conditions  (15)  and  (16)  yield 

n  (a)  =  0  (22) 

n*(0)  »  -tan ( 6  -  j)  .  (23) 

The  functions  n^x')  and  n^x1)  are  obtained  by  integrating  (14)  twice.  The  four 
constant  of  integrations  and  the  value  of  c  have  to  be  evaluated  to  satisfy  the  six 
conditions  (20) -(23).  This  yields  a  system  of  six  equations  with  five  unknowns.  Therefore 
we  cannot  expect  this  system  of  equations  to  have  a  solution  for  any  value  of  Y  other 
than  YQ(B,w). 

The  physical  reason  why  the  problem  does  not  have  a  solution  for  Y  *  YQ  is  that  it 
requires  the  cavitation  number  in  the  sub-bubble  to  be  the  same  as  in  the  main  bubble.  It 
is  to  be  expected  that  the  cavitation  number  within  the  sub-bubble  will  have  some  value 
other  than  Y,  which  we  cannot  prescribe.  Following  the  general  philosophy  of  the  method 
used  by  Vanden-Broeck  and  Keller^  we  shall  introduce  the  unknown  cavitation  number  u  in 
the  interval  c  <  x'  <  a. 
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Integrating  (14)  twice  we  obtain 

ni  (x* )  -  (a2  -  x*  2)  a  (a  +  x'  )log(a  +  x ' ) 

-ja(a  -  x • ) log ( a  -  x 1 )  +  A  +  Bx*  ,  (24) 

n2(x#)  *  (a2  -  x*2)  -  |a(a  +  x’)log(a  +  x') 

•  j  ala  -  x' ) log! a  -  x ' )  +  E  +  Dx  *  •  (25) 

Here  A,  B,  E  and  D  are  the  four  constants  of  Integration.  Using  the  six  conditions 
(20) -(23)  we  obtain  a  system  of  six  algebraic  equations  for  the  six  unknowns  A#  B,  E, 

0,  M  and  c.  This  system  can  easily  be  solved  and  yields  a  unique  solution  for  any  Y  in 

2x 

the  interval  -®  <  Y  <  Yq(B,w).  Typical  profiles  for  8  ■  -j-  are  shown  in  Figures  2  and 
3.  The  value  of  YQ  is  approximately  equal  to  -1.7.  As  Y  decreases  the  size  of  the 
sub-bubble  increases  and  the  size  of  the  main  bubble  decreases.  For  Y  *  W  *  -0.39 
and  the  main  bubble  vanishes.  It  is  interesting  to  note  that  the  present  solution  also 
exists  in  the  interval  YQ  <  Y  <  Y*.  Here  Y*  is  the  value  of  Y  for  which  n^(c)  *  0. 

A  similar  result  was  found  by  Vanden-Broeck  and  Keller'*. 

The  results  are  summarized  in  Figure  4.  The  solution  before  contact  described  in 
Sections  2  and  3  correspond  to  the  interval  YQ  <  Y  <  *•  It  is  represented  by  the  straight 
line  M  -  Y  in  Figure  4.  The  other  curve  in  Figure  4  corresponds  to  the  present 
solution.  It  exists  in  the  interval  -*  <  Y  <  Y*.  Therefore  there  are  two  possible 
solutions  in  the  interval  YQ  <  Y  <  Y*. 
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>.  Solution  vith  an  interval  of  contact 
In  this  section  we  derive  another  solution  for  Y  <  ( 0 * w )  by  requiring  the  bubble 

to  be  collapsed  between  x*  *  f  and  x*  **  g  {see  Figure  5).  We  describe  the  profile  of 
the  bubble  by  the  equations  y*  =  n^(x*)#  0  <  x*  <f  and  y'  »  g  <  x*  <  a.  The 


functions  n.tx')  and  n_(x')  roust  satisfy  the  following  conditions 


I 

I 


i 

I 


Figure  5,  Profile  of  the  bubble  with  one  segment  of  contact  for  Y  *  -3.0  and 
fJ  =  2*/3.  The  vertical  scale  is  the  same  as  the  horizontal  scale. 

The  values  of  f  and  g  are  respectively  0.19  and  0.47. 

The  functions  n^x')  and  n^(x')  are  obtained  by  integrating  (14)  twice.  They  are 
therefore  given  by  the  relations  124)  and  (25).  The  six  constants  A,  B,  E,  D,  g  and  f 
are  found  by  satisfying  the  six  conditions  (26)- (29).  We  note  that  the  present  solution 
can  be  found  with  the  same  cavitation  number  everywhere. 

A  typical  profile  for  0  »  is  shown  in  Figure  5.  As  Y  decreases  the  sizes  of 
the  main  bubble  and  of  the  sub-bubble  decrease.  Furthermore  the  length  of  the  contact 
segment  increases  as  Y  decreases.  For  Y  *  -®r  the  bubble  reduces  to  a  straight 
segment  a  length  a  lying  on  the  x'  axis. 

Finally  let  us  mention  that  the  equilibrium  of  forces  require  the  segment  of  contact 
to  be  a  "film  of  impurities**  characterized  by  a  surface  tension  equal  to  20.  This  is  very 
unlikely  to  occur  in  reality.  Therefore  the  bubble  with  a  segment  of  contact  is  physically 
unrealistic.  However,  this  mathematical  solution  is  physically  relevant  to  describe  the 
deformation  of  an  inflated  membrane.  For  details  see  Vanden-Broeck4. 
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equation  generalizes  the  integrodif ferential  equations  derived  by  Vanden-Broeck 
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and  Keller  .  The  shape  of  the  bubble  is  found  approximately  by  using  the 

2  ' 

slender  body  theory  for  bubbles  presented  by  Vanden-Broeck  and  Keller  .  When 

y  reaches  a  critical  value  Yq(8#<$),  opposite  sides  of  the  bubble  touch  each 

other.  Two  different  families  of  solution  for  y  <  y^  are  obtained.  In  the 

first  family  opposite  sides  touch  at  one  point.  In  the  second  family  contact 

is  allowed  along  a  segment.  The  methods  used  to  calculate  these  two  families 

2  '  4 

are  similar  to  the  ones  used  by  Vanden-Broeck  and  Keller  and  Vanden-Broeck  . 


